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1 Introduction 



Inspired by AdS/CFT correspondence [T], several attempts to construct a holographic 
model of QCD appeared recently, in both bottom-up[2j [3j H] and top-down[5] approaches. 
The morale is to introduce a dynamically generated additional space, which roughly cor- 
responds to the energy scale of the field theory, and try to construct a holographic dual 
model that captures important aspects of the original 4D field theory, such as large N 
expansion, confinement, and chiral symmetry breaking et cetera. The amazing obser- 
vation is that these non-perturbative aspects of strongly coupled field theory, which are 
highly non-trivial to analyze, can be described in simple terms in the 5D holographic dual 
model, sometimes called AdS/QCD mode0. In a sense, the holographic model may be 
viewed as a relevant effective theory for strongly coupled, non-perturbative field theory. 
The counter-intuitive aspect at first sight is that this effective theory necessitates a new 
space dimension generated dynamically. Note that this does not necessarily mean that 
the additional dimension is physical; in view of the 4D field theory, it is simply a derived 
notion that facilitates a magically simple description of important aspects of strongly 
coupled field theory. 

The origin of the simplicity in the 5D dual model is traced back to the large N 
limit as was the case in the original AdS/CFT correspondence in string theory. We 
should thus keep in mind that a usual simple, tree-level calculation in the holographic 
dual model is capturing the leading N contributions, and we are bound to suffer from 
sub-leading correctionsjj Within this limitation, the AdS/QCD model may be taken as a 
new alternative effective theory of low energy QCD. 

The simplest bottom-up construction for AdS/QCD (2J [3], sometimes called the 
hard-wall model, takes a slice of 5D AdS spacetime^l 

ds 2 = — (-dz 2 + dx^dx v r}^ , (1) 

with < z < z m , where z m is the IR cut-off corresponding to the confinement scale, 
which must be fitted to physical observables later. In this paper, we will concern only 
with Np = 2, (u,d) sector of QCD, whose extension to = 3, (u,d,s) sector will be left 
open to a future work. Corresponding to the global chiral symmetry SU(2)l x SU{2)r 
of QCD, we introduce local gauge fields Al and Ar in the 5D AdS slice, whose values 
at z = act as external sources for SU(2)l and SU(2)r currents respectively. More 
precisely, (Ai) 1 ^\ z= q couples to (ji)^ = QlI^Qlj and similarly for SU(2)r sector. The 
chiral symmetry breaking, explicitly by the small current quark mass (Mg^-g^g^. + h.c, 
as well as spontaneously by the chiral condensate E y = (^Rq}^ ^ is holographically 
realized by a VEV of a scalar field X which is bi-fundamental with respect to the gauge 

4 We will later refer AdS/QCD model for phenomenological bottom- up constructions only, but for now 
we mean both top-down and bottom-up. 

5 There is no framework in phenomenological AdS/QCD models yet to do a systematic expansion. 

6 Our metric convention is rj MN = diag(-| ), and we denote 5D coordinates by capital letters 

while 4D coordinates are written with Greek letters. We will follow mostly the convention in Ref.[5J. 
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group SU(2) L x SU(2) R , 



(X(z)) = l -M q z + i£z 3 = {±m q z + \az^ 1 , (2) 

where in the last equality, we assume iso-spin symmetry for simplicity. From the standard 
5D action for the field X, 

S x = J dzj dx A i/G~ 5 Tr (\DX\ 2 + 3|X| 2 ) , (3) 

with DX = dX-iA L X + iXA R , the above VEV breaks SU(2) L x SU(2) R to the diagonal 
iso-spin SU{2)i as expected. The tower of (pseudo) vector and scalar mesons, including 
pions, arise from the normalizable KK modes of the above 5D fields. In addition to the 
parameters appearing in the above, the 5D SU{2) L x SU{2) R gauge coupling constant in 

S g = JdzJ dx i y/G~ s Tr (^--Lf mn F mn ^ , (4) 

specifies the model completely for the two-flavor meson sector. After fitting these param- 
eters, the physical observables are predicted, which compare well with experiments [2l [3j. 

We point out that this phenomenological model in fact encodes essential ingredients 
of the top-down model of Sakai-Sugimoto[5], in that the relevant 5D part of the 10D 
spacetime is approximated by an AdS slice and the world volume gauge fields on two 
branches of D8 branes are identified with Al and A R respectively, whose breaking to the 
diagonal SU(N F )j is achieved by tachyon condensation (X) ^ that joins the two D8 
branes. 

Since the nucleons, i.e. the protons and neutrons, are integral objects in QCD, it is 
pertinent to include them in the AdS/QCD model. An early attempt to see the high spin 
versus mass Regge trajectory was pioneered by Brodsky-Teramond|6J, but recently Ref. [7] 
focused on the lowest spin | nucleons with chiral symmetry breaking effects carefully 
included, and provided us with a concrete AdS/QCD model with holographic nucleons. 
This will be our framework for studying the physics of finite nuclear density. 

We comment that the model in Ref . [7J again mirrors the salient aspects of the top- 
down description of spin | baryons in the Sakai-Sugimoto model [8]. Spin | baryons from 
quantized small instanton-solitons on the D8 brane [HI QUI [H] are effectively described by 
a 5D Dirac spinor field of SU(2)p fundamental. As we identify two branches of D8 brane 
as representing separate SU{2)l and SU{2) R sectors in the phenomenological AdS/QCD 
model (which are combined to SU(2)j at the IR), the 5D Dirac spinor on the two branches 
of D8 brane will corresponds to two separate 5D Dirac spinors in the AdS/QCD model, 
which are fundamental under SU(2)l and SU(2) R respectively. Indeed, this is the field 
content of the model in Ref. [7j. 

Chiral symmetry is believed to be restored at high density and/or temperature, such 
as in the core of neutron stars or in the early stage of fire balls in the relativistic heavy 

7 Iso-spin violating effects from m u ^ md are neglected in this work. 
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ion collisions. Since the discovery of the pulsar [12] , neutron stars have provided extreme 
environments to test the validity of general relativity, especially the gravitational wave 
radiation [13] and the emission of short hard gamma ray bursts [H] from neutron star 
mergers. Many observatories to detect the gravitational waves (such as LIGO, GEO, and 
Virgo) and gamma-ray bursts (such as Swift and HETE II) from neutron star mergers 
are in operation or under development. Neutron stars can also set-up nuclear matter of 
extreme density in which chiral symmetry is restored. Various phase transitions, such as 
kaon condensation [15] . are suggested in the dense core of neutron stars [HI [17]. Even 
though the core of neutron star cannot be observed directly, the maximum mass of neutron 
star which can be observed strongly depends on the equation of state of neutron star [T7j . 
The change in the maximum mass of neutron stars results in the change in the ratio of 
neutron star-black hole binaries to double neutron star binaries, which can be detected 
by the gravitational wave observation [HI [19]. In addition, the inner structure of neutron 
star will be revealed by a detailed pattern of the detected gravitational wave radiation 
from the gravitational wave observatories. Hence theoretical investigation on the physics 
of extreme dense matter is very important in order to understand the physics of the 
neutron stars. 

At high temperature, in addition to the chiral phase transition, QCD deconfinement 
phase transition from hadronic matter to the perturbative (weakly interacting due to 
asymptotic freedom) quark gluon plasma (QGP) is expected. Such QGP states might 
have been formed at the early stages of the evolution of our Universe after Big Bang. 
However, experimental results of the recent relativistic heavy collision (RHIC) indicate 
that the matter after deconfinement (and chiral) phase transition is nothing like a weakly 
interacting perturvative QGP [20l |2T] . Instead, the matter formed after phase transition 
is a strongly interacting quark gluon liquid. These results initiated very active research 
on the physics of hot and dense matter after chiral and deconfinement phase transition 
[221 [231 123 125]. 

Since the properties of nuclear matter cannot be calculated from the fundamen- 
tal QCD Lagrangian due to strong coupling, various effective theories, such as chiral 
perturbation theory, have been introduced to treat nuclear matter phenomenologically. 
Most of the parameters in these approaches were fixed by the experimental values at low 
densities (from vacuum to normal nuclear matter density), at which chiral symmetry is 
spontaneously broken. Hence, the validity of the extrapolation of the effective theories, 
which were built in the symmetry broken phase, to high densities where chiral symmetry 
is restored, may be questioned [TS]. In this respect, recent holographic dual AdS/QCD 
approach which is based on the chiral symmetry at its construction provides a new tool 
to investigate the physics in the symmetry restored phase. 

In this work, we introduce a new way of describing nuclear matter in the framework 
of AdS/QCD model with holographic baryon fields. We turn on the nuclear density in 
mean field approach by a non-zero 5D bi-nucleon condensate, and this subsequently affects 
various other holographic fields, such as X and U(1)b gauge field V, to result in several 
physical consequences that we study in this paper. We point out the relation between our 
approach and the recent approaches by baryon chemical potential[26l [271 1211 [29]. In the 
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spirit of AdS/CFT correspondence, the quark chemical potential in 4D QCD is encoded in 
the boundary value of the 5D U(l) b bulk gauge field, liq^ij) <-> V (x, z) = n q + • • • . This 
is quite similar to a way of introducing various chemical potentials in chiral perturbation 
theory, where the global chiral symmetry is promoted to a local gauge one, and a chemical 
potential is introduced as the time component of the gauge potential [321 E3] . The physics 
behind this approach can be thought of as follows. When we introduce finite density of 
baryons as instanton-solitons, they are energetically pulled down towards IR boundary 
and deposit there. From the Chern-Simons coupling, these IR localized instanton-charge 
density will source the time-like component of the 5D U(l) b gauge field Vo, and the vector 
potential profile of non-zero baryon chemical potential results. A hidden assumption in 
this picture is that the IR-deposited baryons have a 5-function localized wavefunction at 
the IR boundary, and we solve free EOM for Vq in the bulk without describing IR-localized 
baryons explicitly. Although this is correct as a first approximation, our analysis based 
on the spread nucleon profile along z-direction can be viewed as a more precise analysis 
of the situation. We also comment that the interaction terms of Vq an d meson fields, say 
TC, feature the commutator structure, £ int ~ [Vo,H], and therefore the chemical potential 
fi q will be inoperative to neutral mesons, with an exception of Chern-Simons couplings 
studied in [29|. 

After briefly reviewing some details about our AdS/QCD model with holographic 
baryon fields in the next section, we evaluate the density dependence of the chiral con- 
densate and the nucleon mass. We also calculate the mass splitting of charged vector 
mesons in the iso-spin asymmetric nuclear matter. Finally, we re-analyze the in-medium 
nucleon mass in a more phenomenological approach, where we use the density dependence 
of the chiral condensate determined through Hellmann-Feynman theorem and the Gell- 
Mann-Oakes-Renner relation]!]]. In this case, we find the nucleon mass drops drastically 
with increasing density up to the normal nuclear matter number density po(~ 0.16 fm~ 3 ). 



2 The model with holographic baryons 

In this section, we briefly summarize the model in Ref. [7] to set-up the stage. For sim- 
plicity, for N F = 2 meson sector we take the simplest hard- wall AdS/QCD model that 
was explained in the Introduction. To describe spin ^, iso-spin | baryons in this set-up, 
with the nucleons (proton and neutron) as the lowest iso-spin doublet, we introduce 5D 
(Dirac) spinor fields N corresponding to the 4D baryon operators O with the same spin 
and charges in the original QCD. Note that once we introduce a field in the 5D dual 
model, we will actually get a tower of 4D states as we can equivalently describe the sys- 
tem by Kaluza-Klein reduction along our compact extra dimension. The idea behind this 
is that these spin |, iso-spin | 4D operators O in QCD can create/annihilate the whole 
spectrum of 4D baryons B n with the same quantum numbers, 

(B n \O\0) ^0 , n = 0,1,2, ••• , (5) 
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and in the holographic dual model, these -B n 's are packed into a single 5D spinor field 
N (of the same quantum numbers) as the normalizable Kaluza-Klein modes along the 
extra dimension. This is why holography maps a 5D field iVtoa 4D operator O, not to a 
specific 4D state such as nucleons. Therefore, any holographic model is bound to include 
the excited spectrum as well, in addition to the lowest state proton and neutron we are 
interested in. 

Looking at the problem more carefully, one may be puzzled by a few things. In 
the AdS/QCD model side, we have two gauge fields A L and A R for SU(2) L x SU(2) R 
chiral symmetry which are broken to the iso-spin SU(2)i by a Higgs mechanism with 
X. Any field must have a specific representation under SU{2)l x SU{2)r. In QCD side, 
nucleons (and their excitations) form a doublet under SU(2)j and since they are states 
in the broken chiral symmetry phase, it is not at all obvious to even talk about their 
SU (2) L x SU (2) R quantum numbers. 

A guide to way-out may be seen in thinking about how the chiral symmetry breaking 
is realized in the AdS / QCD side. The symmetry SU (2)^ x SU {2) R exists as a definition of 
the theory, and its breaking is achieved by putting a VEV of X; this allows a theoretical 
or fictitious limit of turning off the VEV of X restoring the chiral symmetry. In QCD 
side, this would again be an imaginary limit of chiral symmetry restoration, and every 
operator will have a definite charge under SU(2)l x SU {2) R . This allows us to match the 
quantum numbers of 5D fields to the quantum numbers of 4D operators. We then analyze 
how chiral symmetry breaking affects these 4D operators, and in the 5D AdS/QCD side 
these must correspond to turning on the VEV of the scalar X. Therefore, the effects from 
chiral symmetry breaking in the AdS/QCD model should be realized by possible gauge 
invariant couplings to the field X. 

Although whether or not the above presumptuous chiral symmetry restoration can 
actually be realized is unclear, it may be taken as simply an intermediate tool to identify 
5D field contents and their couplings to the Higgs field X. This chiral symmetry restora- 
tion limit was in fact considered by 't Hooft in his arguing that the massless chiral nucleon 
doublet should exist in this limit to match UV anomalies of SU(2)l x SU {2) R with those 
in IR. More precisely, a left-handed chiral doublet Ol = {pL, n L) T of fundamental rep- 
resentation under SU(2) L , and similarly a right-handed chiral doublet Or = (pR,n R ) T 
of SU(2)r exist in this limit, and the chiral symmetry breaking to the diagonal SU(2)i 
induces a mass coupling OlOr + \\.c. that results in the observed nucleon mass m^v ~ 0.94 
GeV. Hence, we are led to introduce two 5D spinor fields N\ and N2 whose charges under 
SU(2) L x SU(2) R are (|,0) and (0, |), corresponding to the 4D operator Ol and Or 
respectively^ Moreover, we have to ensure that in the chiral symmetry restoration limit 
(X) = 0, we must have a massless left-handed zero mode from N\ and a massless right- 
handed zero mode from N2 in their Kaluza-Klein reduction. They are the massless chiral 
nucleon doublets to match chiral anomalies in the restoration limit. This requirement in 
fact unambiguously fixes the IR boundary conditions for Ni and N 2 . 

8 By Ol,r, we mean the QCD operators ((?l<7l<7l) and {quququ) that create massless chiral nucleon 
states, with a slight abuse of notation compared with the above which means the massless states them- 
selves. 
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The chiral symmetry breaking by (X) 7^ must then introduce a mass coupling 
between N\ and N 2 . The relevant gauge invariant coupling of the lowest dimension is 



C m = -gN x XN 2 + h.c. , (6) 

with a strength g that must be fitted to reproduce the nucleon mass rriN = 0.94 GeV as 
the lowest mass eigenvalue of the KK spectrum. Note that the would-be massless chiral 
zero modes from iVj. and N 2 are lifted to the single massive nucleon state with = 0.94 
GeV, while the mixing between already massive excitations from separate Ni and N 2 will 
split them into a parity doubling pattern of massive excitations. Here, the definition of 4D 
parity involves the exchange of N\ and N 2 , so that the two nearly degenerate states split 
by mixing jVj and N 2 have opposite 4D parity. The parity doubling for excited baryons 
is thus a prediction of the model. 

One last thing to mention is the chirality of the 4D QCD operators Ol and Or and 
its relation to the 5D Dirac spinors N\ and N 2 . One may be puzzled about this since 5D 
spinor doesn't have a chirality. The resolution lies in the fact that in 5D the signature of 
the Dirac mass term flips its sign under parity. The magnitude of the 5D Dirac mass 777.5 
for N\ and N 2 is fixed by an AdS/CFT relation[ 



(m 5 ) 2 = (A-2) 2 , (7) 



with A beingthe scaling dimension of Ol,r, which we take the free theory value A = | 
for simplicityo With a positive 7775 for a 5D spinor N, it can be seen that only the right- 
handed component Nr of N survives kinematically near the AdS boundary z = 0, while 
the left-handed component decays fast enough to become a normalizable modeo The 
non-normalizable Nr near the boundary then couples to a left-handed chiral operator Ol 
in QCD side by NrOl + h.c. Note that this is consistent with the fact that N and Ol 
have the same quantum number, and the possible normalizable chiral zero mode from N 
is left-handed matching to the chirality of Ol- With a negative 7775, the story is simply 
reversed!^] Therefore, the right choice for our model is 777,5 — § for Ni and m 5 = — | for 
N 2 respectively. 

In summary, the AdS/QCD model of spin |, iso-spin | baryons is given by the 



action 



Skin = I dz I dx 4 \/Gc, 



iN^Dm^ + iN 2 T M D M N 2 - ^AyVi + ^N 2 N 2 



S m = J dz J dx A J(h [-gN-yX^ - gN 2 X^N! 



where the covariant derivatives for Ni and N 2 include the gauge group SU(2)l x SU(2)r 
as well as the metric connection, and a single parameter g should be fixed to reproduce 



We believe considering anomalous dimension to A would be an important refinement from the analysis 
in Ref.[7]. 

10 Our convention is i^ipL = +fpL and the mass term is C m = —m^NN. 

11 See also Ref. 31 for the case when both chiral components are normalizable with a specific range of 
TO5. 
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the lowest mass nucleoli mass = 0.94 GeV. By expanding Ni and N 2 in terms of KK 
eigenmodes, it is easy to find the eigenmode equations that must be solved to find the 
mass spectrum of 4D spin ~ baryons. Writing N\(x, z) = fiL(z)Bi(x) + fiR(z)Bu(x) and 
similarly for N 2 (x, z) = f 2 L(z)BL(x) + f2n(z)Bji(x), where -Bl,b are the components of 
the 4D eigenmode spinor B = (Bl, Br) t with mass to be determined, we have 



-m N 





mN T > (9) 



with A (= I in our case) in general. As mentioned before, the existence of the chiral zero 
modes when (X) = requires us the IR boundary condition fiB,(z m ) = f 2 L(z m ) = 0. From 
the meson sector, Ref. [21 [3] found the best fit for (X) = \{m q z + az 3 ) with m q = 2.34 
MeV and a = (311 MeV) 3 , as well as the IR cut-off z m = (330 MeV)- 1 . Then, the only 
remaining parameter of the theory is the dimensionless coupling g, which was found to 
be g = 9.3 in Ref. [7] to reproduce = 0.94 GeV as a lowest mass eigenvalue. 



3 Finite nuclear density in AdS/QCD 

The baryon chemical potential, introduced as a background for the time component of 
the bulk U(1)b gauge field [23, [271 EH EH], is of limited use for some observables due to 
a specific commutator structure of the interaction terms. An example would be its effect 
to the chiral condensate which is encoded in the VEV of the scalar X because DX = 
dX — i[V,X]. In the present work, we introduce the nuclear density through the mean 
field of the 4D nucleon-bilinear, i.e., p s = (ip{x)_iMx)), where ip(x) is the 4D nucleon field, 
and p s is the iso-scalar baryon number density, From the term C m = —gNiXN 2 + h.c, 
we can discuss its effects on the chiral condensate, et cetera. 

3.1 In- medium chiral condensate 

Without nuclear density, the Yukawa coupling C m in Eq. ([H]) is a purely cubic interaction 
term, and it does not enter the equation of motion for the bulk scalar field X or its VEV, 
(X) = Xq. As is well-known, Xq has the following profile 

X (z) = (^m q z + ^a z 3 )l, (10) 

where m q is the current quark mass and o"o ~ (0.330 GeV) 3 is the chiral condensate in 
the vacuum. For simplicity, we take the chiral limit, m q = 0, as our result is expected 

12 Note that the iso-scalar density p s is roughly equivalent to the baryon number density ps = 
(tpT (x)ip(x)) at low density [15] : p s ~ Pb — (-^5-^ V'}: where p is the baryon momentum. 
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not to change much from a small quark mass. As we turn on a finite nuclear matter and 
introduce the non-zero mean field ^N\N 2 + N 2 N\) of the baryon bi-linear, 

(AW + N,N 2 ) = {f 2R - f lR ) =► (f 2R - fl R )p s , (11) 

where N x = fi L {z)i) L {x) + fm(z)?p R (x), N 2 = f 2L {z)i) L {x) + f2R(z)ip R (x), and we used 
the important parity relation 

( hi ) = ( -i o ) ( fT R ) 

which is a consequence of even parity of the nucleons, the interaction (jSJ) modifies the 
equation of motion for X and hence the chiral condensate. We assume that only the 
lowest lying baryons, proton and neutron, have non-zero mean field value of their bi- 
linears, since the higher baryon resonances cost too much energy to form a Fermi sea. At 
very high density, of course, we will have to consider dense matter of higher resonances, 
too. The equation of motion for X is now 

and we solve it numerically with p s ~ p# as a varying parameter. 

It is an important question to specify a physically sensible IR boundary condition 
for X at z = z m , while our chiral condensate a sits in the UV asymptotic X(z) ~ ^cz 3 as 
z — > in the massless limit @ Physically, what drives the chiral condensate in the vacuum 
or the VEV of X would be a strong QCD IR dynamics, which may be attributed to some 
unspecified dynamics on the physical IR brane at z = z m in the present AdS/QCD model. 
We naturally expect that this IR-localized dynamics is not very much affected by the 
presence of our bulk bi-fermion condensate, as we imposed the IR boundary condition on 
our baryons such that their bi-linears vanish at z = z m . Therefore, we are led to identify 
X(z) at z = z m with its value in the vacuum, X(z m ) = |o"o^, as our IR boundary 
condition. We show our numerical result for the ratio of the density dependent chiral 
condensate to that in the vacuum, ^, in fig. CD, where c = Pb/po with p being the normal 
nuclear matter density. We see that the chiral condensate is a decreasing function of the 
nuclear density. 

3.2 In-medium masses of nucleons and charged vector mesons 

In this subsection, we calculate the density dependence of the nucleon mass as well as the 
mass splitting of the charged vector mesons. We first consider the nucleon mass. As seen 
in the previous subsection, the profile of X depends on the iso-scalar density, and plugging 
it back to our 5D mass coupling ([8]), the nucleon wave function equation dHJ) will also get a 
repercussion resulting a shift in its eigenvalue nucleon mass m^, not to mention its wave 

13 This means that the UV boundary condition is — > as z — > 0. 
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Figure 1: Density dependence of the nucleon mass and the chiral condensate normalized with 
respect to those in the vacuum. Here c = pb/po- 

functions. In fact, for a fully consistent analysis, we should solve the equation (fl3|) for 
X and the equation ([9]) for the nucleon wave functions /(i,2)(L,ij) simultaneously, as we 
have a cubic interaction ([H]) between them. For this purpose, we performed an iterative 
analysis until we get a stabilized solution. Therefore our result is numerically reliable, 
and the result is shown also in Fig. [TJ We observe that the nucleon mass decreases in the 
finite nuclear density. 

In light of the results from other conventional QCD effective theories, for example 
Ref. [36, 39, HOj HI], both the in- medium nucleon mass and the chiral condensate have 
been argued to drop by ~ 20% at p B = p (c = 1). Our results for the density dependence 
from Fig. [T] seem a little smaller than these previous estimates. This might be attributed 
to the fact that our analysis is based on various crude assumptions in the AdS/QCD 
model, but it is also worth exploring possible improvements of the model. One thing 
is that we have ignored a back-reaction on the metric due to the finite nuclear density, 
for example, see [Ml ES]- Suppose we incorporated the back-reaction due to the nucleon 
density, then a modified metric would depend on the nuclear density, and would introduce 
additional density dependence into the game other than Fig. [IJ Changing the metric 
would correspond to setting-up a new effective theory around the Fermi surface studied in 
[36] [37] [38] , similar to the spirit of quasi-particles. Instead of delving into this endeavor, 
one possible phenomenological way of incorporating this new density dependence is to 
simply accept the density dependent chiral condensate determined from other low energy 
effective theories, say in Ref.[41j. and to focus on solving the nucleon mass (Q only. This 
will be discussed in the next section. 

We now move on to the mass of vector mesons at finite density. The equation of 
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motion for the time component of the bulk vector field turns out to be 

[f z - ±0,) V*(z) = -gi^(f 1L + f 1R )( Pp - Pn ) , (14) 

where p p {p n ) is the proton (neutron) number density, and the index 3 represents the iso- 
spin quantum number. Here we assume an iso-spin asymmetric environment p p ^ p n , 
which in turn induces an iso-spin chemical potential. This means that once the iso-spin 
density is turned on, Vq will develop the following classical profile, V^ 3 = c\ + c 2 z 2 as 
z — > 0, as a consequence of solving (fl~4l) . where c\ is the iso-spin chemical potential pi 
(= p p — p n ) and c 2 corresponds to the iso-spin density pi{= p p — p n )- Let us denote 
p p = xps and p n — (1 — x)pb- In solving the above equation ffl~4|) . one may take two 
different approaches. One would simply take the density as an input in the UV boundary 
condition, and impose the Neumann boundary condition O z Vq = at IR z = z m without 
breaking the symmetry by boundary condition. This is in accord with our motivation 
that a nuclear density is introduced as a bulk effect, without resorting to an unspecified 
IR-localized baryon dynamics to generate Vq profiles. Another approach would be to treat 
both pi and pi as inputs, as we normally know the one when we are given the other. For 
example, p p>n = \Jm 2 N + (kp n ) 2 and p v ^ n = (kp n ) 2 / (3n 2 ) relate the two by kp n . In this 
case, we start from the given UV z —>■ asymptotic by pj and pi, and solve (fT4l) into the 
IR region. In the present work, we adopt the latter scheme. 

To obtain the in-medium mass of the charged vector mesons, V^, which are linear 
combinations of and V 2 , we consider their equation of motion in the presence of finite 
nuclear density. With a Kaluza-Klein mode expansion, V^(x, z) = g$ J2n fn( z W^ n H x )i ft 
is not difficult to arrive at 

(d 2 z - -d z + mf ± 2F(z)m± + F(z) 2 ^j f±(z) = , (15) 

where F(z) is the solution of the previous equation (|14[) for V^ 3 . Here we impose the fol- 
lowing boundary conditions: /^(0) = 0, d z f^{z m ) = for a normalizable eigen-f unction. 
We focus on the n = 1 mode, which are the charged p-mesons, and the results are shown 
in Fig. [21 Note that we are considering negative iso-spin chemical potential, which is rele- 
vant for the neutron star. As \pi\ increases, the energy cost of putting additional particle 
with a negative iso-spin quantum number into the system will decrease, and the mass 
of the particle will be reduced by \pi\ [53]. Our results in Fig. [2] confirm this physics. 
We remark, however, that as x — > 0, i.e., as \pi\ increases, it is known that the pion 
condensation comes into the story [33J, and the ground-state of the system needs to be 
redefined, and therefore our results in Fig. [2] are no longer valid at small x. 

4 More phenomenological approach 

As mentioned in the previous section, density dependence of chiral condensate in our 
work is rather weak, and this leads to the corresponding weak dependence of the nucleon 
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Figure 2: Mass splitting of charged p- mesons in iso-spin asymmetric matter (x 7^ 0.5). Here x 
is the proton number density fraction, p p = xps- Note that small x corresponds to large 

mass on the nuclear density. If one is simply interested in the nucleon mass in the 
finite nuclear density, it is logically meaningful to simply take the chiral condensate from 
other effective theories that one might trust more, and to just focus on its effects to the 
nucleon mass through the 5D mass coupling ([S]). More explicitly, we take X = C2Z 3 with 
c 2 = cr(ps)/2, where cr(p B ) is the density dependent chiral condensate from the results of 
previous estimates, and we solve ([9]) to re-calculate the in-medium nucleon mass. Another 
possibility may be to use the in-medium p-meson mass, instead of the vacuum p-meson 
mass ~ 770 MeV, to fix the density dependent IR cutoff z m , but we will not pursue this 
in the present work. 

We adopt the model-independent ct(pb) in Ref.jIT], which is valid at low density^ 



and plugging ffTBT) into the KK mode equation for the nucleons, we obtain the in- 
medium nucleon mass as shown in Fig. [31 At normal nuclear matter density, it drops about 
30 %, which seems to agree with the results from other low energy effective theories of 
QCD [36l [39l HQ], HI] . This indicates that the nucleon mass is rather closely related to the 
chiral condensate. 

14 Higher order corrections to the model independent in-medium chiral condensate are studied in 
Ref . [42] . and it is shown that there are at most 15 % corrections up to two times the normal nuclear 
matter density. 




(16) 
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Figure 3: In- medium nucleon mass in the more phenomenological approach. Here Rm = 
M N (p B )/M N {p B = 0). 

5 Summary 

In this paper, we study the physics of finite density of nuclear matter in the framework of 
AdS/QCD with holographic baryon fields. Based on a mean field type approach, we turn 
on the nucleon density by bi-nucleon condensates, and calculate the density dependence 
of the chiral condensate and the nucleon mass. Our result shows that that the chiral 
condensate as well as the nucleon mass drop with an increasing nucleon density. We also 
investigate a mass splitting of the charged p-mesons in an iso-spin asymmetric nuclear 
matter, and the result is compatible with our expectation from a simple physical reasoning. 
We finally study the in-medium nucleon mass in a more phenomenological approach and 
find that the decrease of the nucleon mass in the nuclear matter is strongly correlated to 
the density dependence of the chiral condensate. 

It would be quite interesting to extend our work to the three flavor case including 
strange quarks. Especially, the scaling of the kaon effective mass at high density is very 
important for the neutron star equation of state, because kaon is the least massive boson 
with a strange quantum number. This is left to a future work. 
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